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Abstract. The dependences on z of fractional number densities of H"*" and He"*" ions are 
calculated with a proper allowance for two-photon decays of upper levels of hydrogen and 
parahelium and radiative transfer in intercombination line 2 ^Pi <-> 1 ^Sq of helium. It is 
shown that for hydrogen this gives corrections for a degree of ionization in no more than 
a few percents but for helium this leads to a significant acceleration of recombination 
compared to the recent papers by Seager et al. (1999; 2000) where these effects were 
ignored. 
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INTRODUCTION 

Dynamics of primordial hydrogen and helium recombination was considered in a number 
of papers (see a brief review and references in Seager et al., 2000). However an actuality of 
new more accurate and detailed investigations of this process does not diminish. This is caused 
by a growing accuracy of new measurements of the microwave background radiation (CMBR) 
parameters with the aim to detect contributions of different new fundamental physical factors 
- dark matter, dark energy etc. 

The most detailed calculations of the matter recombination in the Universe were fulfilled 
by Seager et al. (2000) by means of numerical solution of nonstationary equations for the level 
populations of hydrogen (300 levels), Hel (200 levels). Hell (100 levels) and for the number 
densities of electrons, protons, hydrogen negative ions H~ and molecules H2 and H2' jointly 
with an equation for a matter temperature. With all this an average radiation intensity was 
taken to be the Planck function at all frequencies except for the ones in resonant lines for 
which the Sobolev approximation have been used. Collisional transitions were also taken into 
account along with radiative transitions but their contribution turned out to be negligible (as 
was obtained by a number of authors). It should be stressed that for Hel there were taken into 
account singlet states (parahelium) as well as triplet states (ortohelium) . 

The main result by Seager et al. (2000) concerning Hel recombination consists in that it 
goes much more slower than in equilibrium case (according to Saha equation) and slower than 
it was obtained by other authors. And as a main "regulators" of recombination rate appear 
transitions from the second level of parahelium - two-photon ones (2 ^5*0 ^ 1 ^^o) and those of 
in a resonant line 2 ^Pi ^ 1 ^Sq. It turned out that (as in the case of hydrogen) the results of 
multilevel calculations are well described in terms of "effective three level atom" (offered earlier 
for hydrogen by Peebles (1968) and by Zeldovich et al. (1968)) for parahelium with the usage 
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of a corresponding fitting of the total recombination coefficient for the upper levels. In this 
connection Seager et al. (1999) created a symplified code recfast.for to compute recombination 
dynamics in terms of "effective three level model" of hydrogen atoms (levels Is, 2s, 2p + united 
continuum of upper levels) and helium atoms (levels Is^S'o, 2s ^5*0, 2p^Pi + united continuum 
of upper levels of parahelium). 

In the present paper we investigate an influence of additional factors on the recombination 
dynamics of hydrogen and helium (Hel), namely - a contribution of two-photon transitions 
from upper levels of hydrogen and parahelium down to the ground state and radiative transfer 
in intercombination line 2 ^Pi ^ 1 ^Sq of Hel, which were not taken into account explicitely 
in the code recfast.for by Seager et al. (1999) and probably were not taken into account 
in their multilevel calculations (Seager et al., 2000). We have considered how one can take 
into account these factors in a code recfast.for and have obtaied a significant acceleration of 
the Hel recombination. We also have written our own computer code for calculation of Hel 
recombination dynamics with a proper allowance for all factors mentioned above and have 
obtained practically the same results as according to a modified code recfast.for. 



Nonstionary equation for a degree of ionization of a chemical element (hydrogen or helium) 
y = N'^ /N in the uniform expanding Universe can be written in a form 



where is a total number density of atoms and ions of an element, and is a number density 
of its ions, H{z) is the Hubble factor, 2 is a redshift, Ri is a net rate of transitions from a level 
i down to the ground level 1 defining irreversible recombination. To pass from the time scale t 
to the scale of redshifts z we have used an equation dz/dt = —(1 + z)H{z). 

Equation (^J) holds in assumption that the level populations of excited states satisfy to 
the stationary equations of statistical equilibrium with number densities of electrons and ions 
available at the moment. This assumption is really fulfilled because the population of excited 
states is defined by permitted transitions (in the field of blackbody background radiation) 
while the number densities of electrons and ions are defined by much more slower two-photon 
transitions as well as by a "red- shifting" of resonant photons to the longwave region due to the 
Universe expansion. We take into account only radiative transitions because as it commonly 
known (see e.g. Seager et al., 2000) collisional transitions are negligible. Then 



where iVj is a population of a level i, An, Bn and Bu are Einstein coefficients for transitions 
between level i and the first level, Ju is some "average" radiation intensity at the transition 
frequency. 

Up to the present time the two-photon decay of the second level (2 ^5*0 for helium) and the 
"red- shifting" of photons in the main resonace line (2 ^Pi — > 1 ^Sq for helium) were considered 
as the only regulators of the recombination rate (see e.g. the paper by Seager et al. (2000) and 
references therein). However it is clear that similar processes should be taken into account for 
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the upper levels of hydrogen and parahelium as well. Moreover, as it will be shown below, the 
intercombination transition 2 ^Pi — > 1 ^5*0 from the lower state of ortohelium takes a noticeable 
part. 

An expression (j21) for two-photon transitions is rewritten as 
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and it is supposed that the reverse process is a "capture" of two photons of a blackbody 
radiation (with a temperature T) with a total energy equal to a transition energy hvu. Here 
vii is a transition frequency, gi is a statistical wheight of a level i. 

Assuming complete frequency redistribution under scatterings and using a boundary con- 
dition that radiation intensity in a shortwave wing of a line strives to the Planck function we 
have for resonant transitions 
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where (3ik is a probability of a photon escape from a process of scatterings due to an expansion 
of the medium (Universe). It is defined through the Sobolev optical distance Tik- 
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Equations (jlj) and (0) correspond to the Sobolev approximation (Sobolev, 1947; see also a 
review by Grachev (1994), devoted to some generalizations of this approximation), though 
within the framework of available kynematics (a uniform expansion) it gives an exact solution. 
Substitution of eq. (0) in eq. (jH) gives 
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where the wavelength of transition Xu = c/vu. So far as at the begining of He"*" recombination 
(for z ^ 2700) we already have for energies and populations of excited states huu ^ kT 
and Ni <C A^^i respectively then in eqs. (jS)) and (jH)) in the last multiplier one can neglect by an 
exponential term compared to unity and in eq. (jB)) - also by the second term in the denominator 
of the fraction. Then the equations mentioned above are rewritten as 



R?'^ = N,A?,'\l - {N,g,/N,g,)e-'^''^^n 
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where Ni = N — because an overwhelming part of hydrogen an helium neutral atoms are 
in the ground state at the epoch of recombination. 

Thus the problem is reduced to a solution of eq. (Q) with Ri at the righthand side according 
to eqs. ^ and (jH)) while the populations of excited states appearing in the righthand sides of 
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these equations are determined from equations of statistical equilibrium for the current values 
of electrons and ions number densities. It should be noted that such an approach was used 
by us earlier (Grachev and Dubrovich, 1991) under the calculations of hydrogen recombination 
within the framework of 60-level model of atoms. However during an almost all time of helium 
recombination [z = 2700 — 1800) the radiation temperature remains high enough for populations 
of Hel excited states (which are less than 5 eV away from the continuum) to be close to 
equilibrium ones (relative to the continuum) i.e. were defined by Boltzmann - Saha equations 
with the electron temperature equal to the radiation temperature (Tg = T): 

iV, = iV,iV+ ^^_^^-^^^ e"-"/^^% giT,) = {27imkT,f'/h\ (9) 

where Ng is an electron number density, huic is a threshold energy of ionization from a level i. 
Then eqs. (jZ)) and (jH)) are rewritten in the form 

Xfi ri 

where n = {2g+ / gi)giT,){l - y)/{N,y). 

In view of importance of allowance for two-photon transitions from upper levels of helium we 
adduce below a brief derivation of equations for transition probabilities according to Dubrovich 
(1987). The process of simultaneous emission of two photons by excited atoms is well-known 
long ago and is described in text-books (see Berestetskii et al., 1989). Well-known long ago is 
also the role of this process in the hydrogen atom for the continuum radiation formation in the 
interstellar medium (Kipper, 1950) and in the early Universe (Zeldovich et al., 1968). However, 
in these specific cases only one state of hydrogen is taken into account namely - 2s. At the 
same time, as it was shown by Dubrovich (1987), for some values of a medium and background 
radiation parameters similar decays of upper levels can give noticable and in some cases a main 
contribution. Here we will consider this question only in the context of hydrogen and helium 
recombination in the early Universe. 

From quantum-mechanical selection rules it follows that actually we should consider is and 
id states only. An exact expression for a probability of spontaneous two-photon transition is 
written according to Berestetskii et al. (1989) (eqs. (59.28)) in the form 
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where u + u' = uu^is, subscripts il, i'V and Is number initial, intermediate and final atom 
states {i is a main quantum number, / is an orbital moment), {dQ)ii^i'v is a matrix element of a 
dipole moment, a and (5 number spatial components of a dipole moment vector, u and v' are the 
frequencies of emitted photons, h is the Planck constant, c is the velocity of light. Sharp maxima 
in eq. (jl2j) under v 01 v' = Vai correspond to resonances of a cascade transition from an excited 
level down with an emission of photons of a discrete spectrum of an atom. Corresponding 
to them very large transition probability leads in specific conditions of quasi-equilibrium with 
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a blackbody radiation to a very high probabihty of a reverse capture of the same photons. 
Strictly speaking this question must be learned more thoroughly because "not entirely" resonant 
transition, but not far from a resonance, can give contribution under the scheme of "a Lyman 
quantum escape into the wing". In principle this will lead to an additional speeding-up of 
recombination. However a thorough analysis requires a radiative transfer solution which we 
intend to obtain in the next paper. Here we will consider only transitions giving continuous 
distribution of emitted photons i.e. we shall assume that un — u' ^ u. So an obtained rate of 
recombination can be regarded as a lower estimate. 

Because we are interested in a two-photon transition to a final state with a zero orbital 
moment (s state) then / may be equal only (s state) or 2 {d state) according to selection rules 
for dipole transitions. For both cases the value of /' may be equal only 1 {p state). In this case 
an expression for W can be simplified significantly. It can be simplified still more if to notice 
that according to summation rule for dipole transitions [Berestetskii et al., 1989, eq. (52.8)] 
nearly 90% of contribution is due to transitions with i' = i [Berestetskii et al., 1989, eq. (52.6)]. 
As a result we come to an expression for a matrix element structure well-known for the 2s — Is 
transition in a hydrogen atom. In our case to calculate W we need only in a proper allowance 
for the frequency difference i.e. we must multiply A2s,is by and sum up the two ways 

of decay (from s and d sublevels) with a proper allowance for their statistical wheigts. Finally 
we obtain for hydrogen the following expression: 

Wn = grA^,''^ = 54 • A^sm ■ [(^ - l)/(^ + iTilU' - 41)^. (13) 

For large i we have approximately Wh = 89^ s^^. The growth of Wu with the level number 
takes place actually up to a some value i. This is caused by an existence of the limit of the 
dipole approximation applicability. Namely the wavelength of emitted photon (which is in our 
case ~ 2/i/ji const) must be greater than the size of an excited state orbit (~ i^). For 
hydrogen i^^x ~ 30 (Beigman and Syrkin, 1983). 

Similar consideration can be carry out also for decays of hydrogen-like states of Hel. For 
i > 6 — 7 such an approximation is quite rightfull for the matrix elements i{s, d) — ip. It 
becomes sufficiently accurate for our aims also for the square of the matrix element ip ^ Is 
if we introduce correction multiplier 1.15 - 1.20 which follows from comparing of oscillator 
strengths of these transitions for Hel and for hydrogen. And of course eq. ()13|) must be 
renormalized again with a proper allowance for the frequency of emitted photons i.e. W for 
hydrogen must be multiplied by (24.6/13.6)^ = 19.4. Finally we have for W-aei- 

Wuei = 9^Af^^ = 1045 -A-iii- l)/(z + l)f\lli' - 41) A. (14) 

A limiting condition for i is here the same as for the case of hydrogen. If dipole approximation 
is applicable (for i < 40) then one should take A = 10 in eq. (jl4|) . Otherwise there is an 
uncertainty connected with poorly known (both theoretically and experimentally) dependence 
of A on a level number. Most probably a contribution of these levels {i > 40) is not very 
large. Approximately it can be taken into account taking A = 12 s~^. However for precise 
measurements of the power spectrum with the aim to obtain information about weak but very 
important factors an additional investigation of A is necessary. 

An influence of two-photon decays of upper levels on the dynamics of recombination must be 
much more significant for Hel than for hydrogen. This is caused by two circumstances: first - 
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by much more absolute value of W and, second, - by the fact that relation between populations 
of 2s-level and Rydberg levels significantly differs since a ratio of these populations contains 
the Boltzmann factor exp{—hi'i2/kTr). For hydrogen this factor is equal approximately 3 • 10~^ 
while for Hel it is larger approximately at 85 times since an absolute values of energy differences 
for Hel are approximately the same but the temperature at which it recombines is significantly 
higher. A contribution to a rate of destruction of "superfluous" Lyman quantums is defined by 
a product of a population on a decay probability. 

We wrote computer programme to solve eq. for helium with the terms in the righthand 
side of the form ()10|) and (fTT|) and along with the two-photon transitions from the second 
level of parahelium (i = 2s ^ 1) we also took into account transitions from upper levels 
{i = 6 — 40) which were considered as hydrogen- like (Rydberg ones). For their energy counted 
from a threshold of ionization we take huic ~ IRy/z^ and for Einstein coefficients we use 
eq. (I14j) . Moreover along with the resonant transition from the second level of parahelium 
i = 2p^Pi <-> 15^5*0 (thereafter we wright for the sake of simplicity i = 2p ^ 1) it was also 
taken into account spin-forbidden one-photon transition from the second level of ortohelium: 
2p^Pi ^ Is ^5*0 for which we use the value of Einstein coefficient ^2 3pi,ii5o = 233 s~^ according 
to Lin et al., 1977. It should be stressed that for the degree of ionization of hydrogen which 
is contained in equation N^, = + N^^ we have used Saha equation (equilibrium ionization) 
which holds for sufficiently large z where the main Hel recombination takes place. 

We also introduced corresponding additions into the programme recfast.for by Seager et 
al. (1999) in which two-photon transitions were taken into account only from the second level 
[i = 2s ^ 1) and Einstein coefficients ^2^"^] for hydrogen and parahelium (Ah and Ahbi in 
notations by Seager et al., 1999) Ah = 8.22 s~^ and Anei = 51.3 s~^. To take into account 
two-photon transitions from upper levels it needs evidently according to eq. (fTnjl to make a 
replacement 

A ^ A + ^ ^e'^^''--"^^'^)/*^^. (15) 

i=io 

Further, in the programme recfast.for a contribution of transition i = 2p ^ 1 for hydrogen and 
parahelium is described by the factors Ku and K^ei respectively where K = Xi^2p/[8'^H{z)] so 
that a contribution of other lines (for one-photon transitions) can be taken into account (in 
accordance with eq. ()11|) ) by the following replacement: 
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and for the main resonant transition during the whole recombination an optical depth ri 2p ^ 1 
so that the corresponding exponential term can be omitted. 

As the computations have shown the results obtained under our own programme and under 
the programme recfast.for modified as pointed above are practically coincide (for the values of 
parameters of the Universe model accepted now). 

RESULTS OF COMPUTATIONS 

Parameters in the problem are equilibrium temperature of the microwave background Tq, 
Hubble factor Hq, the ratio of the total density to the critical one fitotah fhe rati of the baryon 
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Figure 1: The profiles of fielium ionization degree. Numbers near tfie curves: 1 - only "main" 
transitions 2s ^ 1 and 2p 1 of paralielium, 2 - main + transition in a line of ortohelium, 3 - 
main + transitions from upper levels of paralielium [A = 10 s~^), 4 - main + both additional 
ways (continuous line - A = 10 s~^, broken line - A = 12 s^^), 5 - equilibrium case (Saha 
equation). Values of parameters: fitotai = 1, = 2.728 K, Hq = 70 (km/s)/Mpc, = 0.04, 
^DM = 0.26, = 0.7. Helium content (by mass) Y = 0.24. 

density to the critical one ^b, the ratio of the dark matter to the critical one Qbm, the ratio 
of the density caused by the A-term to the critical one Q\ at the present day epoch and the 
content of the primordial helium (by mass) Y. Test computations were fulfilled for Tq = 2.728 
K, fitotai = 1,Y = 0.24, Ho = 70 (km/s)/Mpc, = 0.04, r^oM = 0.26 and = 0.7. 

Fig. 1 shows the results of computations. Numbers 1, 2, 3 and 4 correspond to a succesive 
allowance for different ways of irreversible helium recombination: 1 - only "main" transitions 
2s 1 and 2p 1 of parahelium, 2 - main + transition in ortohelium line, 3 - main + 
transitions from the upper levels {i = 6 — 40) of parahelium for A = 10 in eq. (fT^ . 4 - 
main + both additional ways (continuous line - for A = 10 s~^ in eq. ()14|) . broken line - for 
A = 12 s~^, 5 - equilibrium case (Saha equation). One can see from the Fig. 1 that succesive 
allowance for additional ways of irreversible recombination significantly enhances the rate of 
Hel recombination. At the same time our results practically coincide (undistinguished on the 
Figure) with the ones obtained under the programme recfast.for modified by the way stated in 
the preceeding section. 

With the aid of the programme recfast.for we took into account additional ways of irre- 
versible recombination for hydrogen too. These ways are two-photon transitions from the levels 
i = 3 — 40 down to the first level and radiative transfer in the corresponding Lyman lines. The 
first way was taken into account according to eqs. (fT3j) and (fT^ and the second - according to 
eq. (fTBj) while exponents in round brackets in this equation were omitted because the Sobolev 
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Figure 2: An influence of additional transitions on the profile of hydrogen ionization degree 
yu = N^/Nii. On the ordinate axis is = 2{y^'^ — y^'")/{y^'^ + y^""), where y^^ corresponds 
to "main" transitions 2s 1 and 2p ^ 1 only, and |/Jf™ corresponds to main + transitions 
from upper levels. Values of parameters are the same as for the Fig. 1. 

thicknesses in Lyman lines Tu ^ 1. It turned out that the first way leads to decrease of hydro- 
gen ionization degree by no more than 4.2% and the second - by no more than 1% and in sum 
- by no more than 4.3% (for the same values of the model parameters as in the case of helium). 
Results are in the Fig. 2. It should be stressed that we have fulfilled in addition similar com- 
putations using our own programme (Grachev and Dubrovich, 1991) based on solution of eq. 
(jT)) jointly with equations of statistical equilibrium for the level populations of hydrogen atoms 
within the framework of 60-level model of an atom and we have obtained practically the same 
result as in the Fig. 2. 



Computations of hydrogen and Hel recombination dynamics are made with a proper al- 
lowance for two-photon decays of upper levels of hydrogen and parahelium and radiative trans- 
fer in intercombination line 2 ^Pi 1 ^5*0 of helium. It is shown that this leads to changes of 
hydrogen recombination rate accessible to discover under the programme "PLANCK". Ob- 
tained results are important for a correct evaluation of small factors defined by an existence of 
dark matter, baryon and nonbaryon parts of the substance mass in the Universe, relativistic 
(at the moment of recombination) particles - neutrino and possibly axions and other low-mass 
low-interacting particles. It is shown that an allowance for new factors of destruction of suprae- 
quilibrium Lyman photons significantly accelerates Hel recombination compared to predictions 
according to the papers by Seager et al. (1999; 2000) where these effects were not taken into 
account. This is also important for the correct estimation of the helium role in the hydrogen 
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recombination. Moreover it is very important for determination of profiles and intensities of hy- 
drogen and Hel recombination lines. In particular an important role of high levels in irreversible 
recombination leads to an appearance of absorption lines in the CMBR spectrum caused by 
transitions in Balmer lines of hydrogen and in different more complex lines of helium. This 
moment is principally important from the point of view of detection and identification of such 
lines. As a next step in this direction it is necessary to consider a role of weakly nonrcsonant 
transitions from upper levels which can lead to an additional acceleration of recombination for 
hydrogen as well as for helium. 

As concerns an influence of our refinements of Hel recombination dynamics on the theoretical 
power spectra of microwave background radiation it can be estimated very roughly from Fig. 17 
in the paper by Scager et al. (2000) where their results are compared with the results by Hu et 
al. (1995) obtained actually assuming equilibrium Hel ionization (according to Saha equation). 
So, since our curve of variation of Hel ionization degree lies practically "in the middle" between 
equilibrium one (according to Saha equation) and that from Seager et al. (see Fig. 1) then one 
can expect that with our refinements the discrepancy with the results by Hu et al. will become 
2 times smaller (e.g. for multipole / = 1500 it will be 1% instead of 2%). At the same time an 
uncertainty of A (see Fig. 1, curves 4) estimated approximatly at 20% should give the same 
relative uncertainty for correction of theoretical power spectrum. 
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of Russian Federation for supporting of Leading scientific schools (grant NSh-1088.2003.2), 
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